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Abstract 

We attempt to make a direct measurement of the weak lensing signal from the WMAP 7-year 
data. We apply the real-space implementation of the optimal quadratic estimator on the maps 
produced by the W-band Differencing Assemblies. We obtain a weak lensing amplitude parameter 
of Ai = 0.98 ± 2.50 after correcting for several sources of bias. The significance of the estimation 
improves to A^ = 1.55 ± 0.93 when information from WMAP simulations is used. Despite the 
demonstrated insensitivity of the real-space estimator to uncorrelated noise, we conclude that this 
detection is not statistically significant. We expect that a full-sky higher-sensitivity experiment 
such as Planck will allow us to make a more significant measurement. 
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I. INTRODUCTION 



Gravitational lensing of the cosmic microwave background (CMB) radiation is a source 
of information on the interaction between the CMB photons and the fluctuations in the 
matter density distributed along their trajectory. This interaction causes the deflection of 
the CMB photons by the gradient of the lensing potential if), which results in the creation 
of correlations among initially uncorrelated modes. Since the matter fluctuations are in the 
linear regime on large scales, the deflections are small and consequently the total deflection 
can be described as a sum of multiple deflections along the unperturbed trajectory (Born 
approximation). The correlations appear in the form of non-diagonal couplings between 
gradients of the temperature and gradients of the lensing potential, hence as a non-linear 
signal which becomes important at small scales. By measuring these correlations, we can 
estimate the lensing deflection cx = Vif> and consequently the lensing convergence k, = 
—V 2 ip/2 which is a direct measure of the projected matter density. 

Detections of a lensing signal in the CMB are usually summarized in the lensing amplitude 
parameter A L which describes the amplitude of the detected lensing power spectrum relative 
to the theoretical prediction. Thus, A L = indicates no lensing, while A L = 1 indicates the 
standard ACDM model and General Relativity lensing theory. Recently there were reported 
direct detections of the lensing signal from experiments designed to look for non-linear signals 
in the CMB, yielding = 1.16 ±0.29 from the Atacama Cosmology Telescope (ACT) data 
[1], and A L = 0.90 ± 0.19 from the South Pole Telescope (SPT) data [2J. Although the 
Wilkinson Microwave Anisotropy Probe (WMAP) was not designed to look for non-linear 
signals, a detection from the WMAP full-sky data would provide a full-sky measurement of 
the lensing signal and complement the other detections. 

We anticipate some difficulties in this detection. Since the experimental noise in the 
WMAP data overwhelms the lensing signal, a statistically significant measurement should 
not be possible and indeed the first detections were made indirectly by correlating the 
results with galaxy counts [31 0]. There is however a reported direct detection of the lensing 
signal from the WMAP data using a kurtosis estimator, finding Al = {0.96 ± 0.60, 1.06 ± 
0.69,0.97 ± 0.47} from the data of the V, W and combined V+W bands respectively [5], 
whereas a more recent less-significant direct detection from WMAP data, using the optimal 
quadratic estimator [UJ, found A L = 1.27 ± 0.98 [7j. 

Here we base our measurement on an estimator with demonstrated insensitivity to uncor- 
related noise [E] and good handling of excision of points [9], which makes it a reliable tool for 
the extraction of the lensing signal from the WMAP data. Similarly to the optimal estima- 
tor, this estimator combines two filters of the temperature map weighted by the inverse of 
the variance. Instead of acting in Fourier space and using all the available sky for the recon- 
struction of the lensing potential at any point, this estimator acts in real space and uses the 
neighbouring points only, which presupposes a truncation of the kernel that convolves the 
two filters. Although truncating the kernel renders the estimator slightly suboptimal, it has 
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a strong physical motivation. Since weak lensing is a small-scale effect coherent over scales 
of a few degrees [10j, information on the lensing potential at a given point is contained in the 
points within the degree-scale, so we expect little gain in considering the information over 
the entire sky [TT]. Thus a direct detection of a lensing signal from the WMAP data is both 
a test of the robustness of the real-space estimator and a comparison of the performance 
with other estimators. 

The manuscript is organized as follows. In Sec. [IT] we introduce the weak lensing estima- 
tor. In Sec. Ill we justify our choice of data and describe the simulations. We then proceed 
to measure the weak lensing signal. In Sec. IV we use synthesised map so as to simulate a 
control experiment, test the bias removal and setup the calibration. In Sec.[V]we use patches 
from cartesian projections of the full-sky maps. The results are summarized in Sec. |VI| 



II. ESTIMATOR OF THE CMB WEAK LENSING 

The deflection of the CMB photons by the gradient of the lensing potential amounts to 
remapping the unlensed temperature anisotropies T to the lensed T in the direction by 
T(0) = T{0 + V^), where the lensing potential ip is the projection on the image plane 
of the gravitational potential of the large-scale structure. For small deflections, the lensed 
temperature power spectrum is 

(f(£')f(£ - £')) = (2n) 2 5(l)C e , + (2tt) 2 [£ ■ £'C e +£■(£- £')C\^\] 1>(£) + 0(^ 2 ), (1) 

where Cg is the unlensed temperature power spectrum. Both the unlensed temperature 
anisotropies and the large-scale structure are assumed to be Gaussian random fields. The 
remapping of one Gaussian field by another induces non-Gaussianities. The lensing signal is 
encoded in the off-diagonal correlations as a coupling between the photons and the gradient 
of the lensing potential. An estimator that captures the information in the coupling consists 
of a quadratic combination of the temperature map optimized so as to yield the minimum 
variance [6j. This estimator will also pick up the diagonal correlations which exist in the 
absence of lensing and which we must remove both from the estimated map and from the 
estimated power spectrum. The diagonal correlations have a much larger amplitude than 
the lensing signal and are estimated from simulated unlensed temperature maps. 

We measure the isotropic component of the convergence tensor, defined as k = — {d 2 + 
dy)ip/2. We use the real-space estimator introduced in Ref. [H] which is based on the convo- 
lution of the square of the lensed temperature map T with a local kernel W as follows 

k(0) = J d 2 + J d 2 0_f(0 + + 0_)f(0 + -0_)W(0,0 + ,0_). (2) 

The real-space kernel W(0, + , 0_) is given by 

W(0,0 + ,O„) = J ^exp[i£ + -0] J ^ exp {- z (£ + .0 ++ £_.0_)]W(£ + ,£_) (3) 
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where the corresponding Fourier-space kernel W (£+,£-) is related to the minimum- variance 
(to leading order) weight function Q(£,£') by 

W(£ + ,£-) = Q(£,£') (4) 

upon the coordinate transformation £ = £ + , £' = (£ + +£_)/2. (See Ref. [8] for details). Here 



1 £■ £'C £ > +£■(£- £')C ]e i_ e] 



Q(£,£') = -Af e --^ — -z-pr- ' (5) 

4 2 [Q, + iV>] [C| £ _£'| + A^|£_£'|] 

where Ni is the noise power spectrum and A/_ is the variance of the estimator defined by 

= /" i [i-^+i-(i-Oq*-«i] 2 

£ 7 (2tt)2 2 [C t , + N £/ ] [C\ t _g\ + N l£ _ £/{ ] ' 

In the ^-integration we sum modes 2 < £, £' < 1200. The validity range of the real-space 
estimator is defined by the largest and smallest lensing modes that can be probed by the 
kernel. Hence, the lower £-limit is determined by the size of the maps, whereas the upper 
£-limit is determined by the size of the kernel. 



III. DATA AND SIMULATIONS 



The relevant information for the reconstruction of the lensing potential lies on the angular 
scales close to the resolution scale of the map, which is set by the beam size. Hence we opt to 
use the maps produced by the four W-band differential arrays (DA) for having the smallest 
beam size 1 . To minimize the contamination of the lensing signal by spurious sources of non- 
Gaussianity, we use the coadded, foreground-removed temperature maps [TJ] at resolution 
r9 2 . Moreover, to exclude foreground-contaminated portions of the sky at the location of 
the galactic plane and point sources, we use the pixel mask KQ85 which keeps 78% of the 
sky area 3 . 

We conceive different experiments based on the specifications of the four DAs of the 
WMAP W-band 4 . For each differential array, we conceive two experiments: a) the exper- 
iment denoted by "theo," which consists of lensed CMB maps synthesised in the flat-sky 
approximation using the temperature (TT) and lensing deflection (aa) power spectra com- 
puted with CAMB [13] (on scales £ < 1200, for a fiducial flat ACDM cosmology consisting of 
h = 0.7, tt cdm h 2 = 0.122, tt h h 2 = 0.0226, n s = 0.96, a 8 = 0.85, r rci = 0.09 and T = 2.725 K), 
and b) the experiment denoted by "obs," which consists of temperature maps derived from 
HEALPix 5 [14J cartesian projections of the full-sky map produced by WMAP. The specifi- 
cations of the experiments are summarized in Table [T| 

1 http:/ /lambda. gsfc.nasa.gov/product/map/dr4/m_products.cfm 

2 http://lambda.gsfc.nasa.gov/product/map/dr4/maps_da_r9_iqu_7yr_get.cfm 

3 http:/ /lambda. gsfc.nasa.gov/product/map/dr4/ptsrc_catalog_get.cfm 

4 http:/ /lambda. gsfc.nasa.gov/product/map/dr4/beam_xfer_get.cfm 

5 http://healpix.jpl.nasa.gov/ 
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Table I: The specifications of the experiments used in this study. The experiments denoted 
by WMAP^° use the power spectra computed with CAMB and detector noise characteristic of 
the Wi DA. The experiments denoted by WMAP^ S use cartesian projections of the corresponding 
Wi DA measured by WMAP. 

IV. MEASUREMENT OF THE CONVERGENCE FROM THE CMB POWER 
SPECTRUM 

A. The input maps 

We first compute the convergence map from simulations of the four DAs of the WMAP W- 
band. We synthesise lensed CMB temperature maps using the theoretical power spectra (TT 
and eta) computed with CAMB and combining the beam window function and detector noise 
characteristic of each DA of the WMAP W-band. (See Appendix B of Ref. [8] for details.) 
We compute the noise power spectrum from a full-sky map synthesised as White Noise x 
Co/v^Vobs, where cq is the rms of the noise per number of observations and iV b s is the 
corresponding map containing the number of observations per pixel. From the noise power 
spectrum, we synthesise a noise map which we add to the beam-convolved, lensed CMB 
map. We call these experiments WMAP^°. 

For each experiment, we also generate unlensed CMB temperature maps from the theo- 
retical TT power spectrum computed with CAMB. The unlensed CMB maps are generated 
using different realizations of the CMB map for the same cosmology, thus simulating dif- 
ferent realizations of the temperature and noise maps. The unlensed CMB maps serve to 
model the Gaussian contribution of the CMB to the four-point correlations measured by 
the estimator of the lensing power spectrum which is the major contaminant of the lensing 
signal [15] . 

B. The convergence maps 

The convergence map is reconstructed by convolving in real space the minimum variance 
kernel W(9, 9 + , 9 J) with the square of the temperature map from which we want to estimate 
the weak lensing signal. Both the lensed CMB power spectrum and a fiducial power spectrum 
for the unlensed CMB power spectrum enter in the calculation of the kernel. We generate 
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Figure 1: Kernel in real space for the four DA's of the WMAP W-band (theo). We plot 
the dominant {m = 0) eigenfunction only. The levels denote the orders of magnitude as fractions of 
the maximum value, delimited by {1, 10 _1 , 1CP 2 , 10 -3 , 10 -4 }. The kernels in all panels are derived 
from the theoretical power spectra combined with the beam window function of each DA. Top left: 
Wi, Top right: W 2 , Bottom left: W 3 , Bottom right: W 4 . 

the kernel for each DA. In the calculation of the kernel we use the lensed power spectrum 
computed with CAMB combined with the window beam function and the detector noise 
characteristic of each DA. 

The kernel can be expressed as a series of eigenfunctions, the dominant being the m = 
mode which we plot in Fig. [I] (See Ref. [5] for the calculation of the kernel in real space.) 
The plots show that the kernels for the four experiments have the same size and similar 
structures within the 1% level, differing only in the smaller-scale structures (of order < 0.1°) 
at sub-percent level, so we expect that the corresponding lensing reconstructions have similar 
ranges of validity. In all cases, for the numerical implementation we truncate the kernel at 
^kernel = 1-1°, which is the spatial extent which corresponds to the 1% level. As previously 
tested, keeping sub 1% levels does not improve the reconstruction of the power spectrum 

The convergence map k reconstructed from a lensed CMB map contains both the lensing 
signal km and a spurious contribution /tu=o due to the unlensed CMB map. The zero-point 
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function of this contribution amounts to an offset which we remove by subtracting the pixel 
average of the reconstructed convergence map k. The resulting convergence map k — (&) pixcl 
is unbiased in the sense that its expectation value is equal to that of Kup. (See Refs. [SJ |H] 
for details.) 



C. The convergence power spectrum 



After correcting for the bias of the reconstructed convergence map, spurious Gaussian 
correlations due to the unlensed CM B still persist which add to the non-Gaussian correlations 
due to the lensing signal. The power spectrum of the convergence due to the lensing signal 
is obtained after subtracting the Gaussian contribution off of the power spectrum of the 
reconstructed convergence map. To estimate the Gaussian contribution, we take the average 
of the power spectrum of the convergence maps reconstructed from various maps of the 
unlensed CMB. The unlensed CMB maps are synthesised from the same theoretical TT 
power spectrum while varying the seed. We estimate the unbiased power spectrum of the 
convergence C^* by considering the weighted average of C^^° over the different CMB 
realizations 



num maps 



CMB 



a 



K|i/>=0 



num maps 



w i I 2^ 



W; 



|Vl=0 



(7) 



and subtracting it off of C#. The weights are the inverse of the variance of the 
estimated power spectrum 

1 r / , N 12 



sky 



Var[£] + [C t 



*|t/>=0 



(8) 



and Var[£] is the variance of the estimator given by Eqn. (|6]). The unbiased power spectrum 
estimated by this method is thus given by 



' CM B ' (9) 

which is a difference between two very similar quantities. The total variance is given by 



Var [C^] = Var [C£] + Var [C^ =0 ] . 



Here 



Var[Cf] = Var est [C;]. 



Var[C* w=0 ] 



num maps 



J2 (l/Var est [C; w=0 ] 



Var s [C* w=0 ] : 



(10) 

(11) 
(12) 



where Var est \Cf\ is given by the analogous to Eqn. (Jsl) upon the replacement k\q, =a — > k^ 



and Eqn. (|12|) includes the sample variance about the average value of °, Var s [C^ °1 . 
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The estimate of the convergence power spectrum averaged over the four experiments 
WMAP{f is plotted in Fig. g (left panel). The sample variance becomes subdominant 
in the entire £-range with the increase of the number of maps. The detection error is 
dominated by the estimator noise set by the first term in Var est , which depends on the 
experimental resolution and detector noise. We observe, however, that the result contains 
a large excess of power as compared to the expected signal, the input convergence power 
spectrum also shown in Fig. [2j In contrast, for simulations of the higher signal-to-noise 
Planck experiment, the quantity {C^ =0 ) CMB was found to be a good estimator of the 
Gaussian bias as demonstrated in Refs. [8j |9]. 

A source of additional bias are couplings of lensing modes which add spurious correlations 
to the estimated power spectrum from the four-point correlation function of the lensed CMB 
[15]. We expect to get a much better agreement if we are able to reduce the uncertainty in 
the CMB map since the couplings add incoherently [15] . This was not required for the signal- 
to-noise ratio of Planck, but it is needed for WMAP. We then produce various lensed CMB 
maps, all from the same CAMB lensing and temperature power spectra, and again generate 
unlensed CMB maps using the same realizations of the lensed CMB maps. Throughout the 
different realizations we generate different realizations of the noise map. 

We estimate the unbiased power spectrum of the convergence by taking the difference 
of the weighted averages over the different CMB realizations 

~ syt 1 cmb ~ syt /CMB- y 16 ) 

Here 

num maps num maps 

( c ")cmb= E WW/ E < ( 14 ) 

i i 

with the weights wf given by the inverse of the variance of the power spectrum Var^C^]^, 



and (C^ °) CMB given by Eqn. (7). The total variance of C™ 1 * is given by Eqn. (10), where 
Var \Cf\ is given by 



Var[C;] 



num maps 



J2 (l/Var^J 



Var s [C;] (15) 



and Var[c£ |l/ ' °] is given by Eqn. (12). 
When applying the estimator o 



Eqn. (13), we indeed obtain a convergence power spec- 



trum, also shown in Fig. [2] (left panel), that closely reproduces the input power spectrum 
within the validity range of the real-space estimator. 

As a measure of the statistical significance of the detection, we compute the best-fit 
lensing amplitude parameter Al- We fit the estimated power spectrum with functions of 
the form Al C™, where C™ is the input fiducial convergence power spectrum. The x 1 - 
minimization yields an analytical expression for Al, which is essentially an inverse-variance 
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Figure 2: Estimated convergence power spectrum from CMB simulations of the 
WMAP's W-band DAs. Left panel: The gray diamonds are the convergence power spec- 
trum computed as the difference between one lensed CMB map and the inverse-variance average 
over different realizations of unlensed CMB maps. The black squares are the convergence power 
spectrum computed as the difference between the inverse-variance average over different realiza- 
tions of paired lensed CMB and unlensed CMB maps. For each DA, the maps are synthesised from 
the theoretical power spectra combined with the beam window function and noise power spectrum 
of the corresponding DA. The dotted line is the theoretical convergence power spectrum computed 
with CAMB. The dash-dotted line is the inverse-sum of the variance of the estimator Var[Ki^] for 
each DA. Right panel: The black line shows the ratio between the two quantities in the left panel. 
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weighting of the ratio Cg ut /C™ summed over all scales i [16J, 



E,Q n Q° u VVar[Cr] 



£,(C7) 2 /Var[C; 



out! 



(16) 



The error of Al is computed by propagating the error of Eqn. (16), yielding 



£ £ (C7)7Var[C< 



out] 



(17) 



We compute Al for each DA, (A L )w i} for the average of the four DA's auto-correlations, 



(A? c 



)(WiW, 



}, and for the full combination of DA's considering all auto and cross-correlations, 



/ /ltheo\ 

(Al ){WiWj)- 

The results for the four experiments, separately and combined, in the validity range of 
the estimator, are presented in Table |TT} For the combined result we find Al = 1.55 ± 0.93, 
which is similar to the findings in Ref. [7j. 
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(Af™) Wl (Af™) W2 « eo )^ 3 « e °)m 


/ /itheo\ ( /itheoN 

\ A L )(WiWi) \ A L )(WiWj) 


t E [6, 262] 


1.49 ±2.95 1.49 ±3.73 1.55 ±3.90 1.59 ±3.45 


1.53 ± 1.72 1.55 ±0.93 



Table II: A^co f rom eac h W-band DA and from the combined W-band DAs. 



V. MEASUREMENT OF THE CONVERGENCE FROM THE CMB MAPS 



We now turn to the problem of detecting a lensing signal in the WMAP 7-year data. 
Having one full-sky realization only of the lensed map, we cannot apply the estimator defined 
in Eqn. (13). We can however use the simulation results to infer an empirical calibration. 
We define the empirical calibration as the ratio between the two estimated in Eqns. ([9 
and (13), and plot it in Fig. |2] (right panel). 



A. The input maps 

We then proceed to compute the convergence map from the temperature maps produced 
by the four DA's of WMAP W-band. We use the foreground-reduced maps and apply the 
KQ85 mask which keeps 78% of the sky area. We then make cartesian projections of the 
full-sky masked maps into patches of the same size (56° x 56°) as the simulated maps in 
the previous section. The projections cover the whole sphere in steps of 45 degrees along 
the latitude and along the longitude. In a similar way, we also project the mask and the 
map containing the number of observations per pixel for each DA which we use to generate 
the noise patches. We call these experiments WMAP^. In Ref. [17], a similar strategy is 
adopted to estimate the bispectrum from a full-sky map of a Planck simulation. 

For each patch we generate unlensed CMB maps. These maps consist of the same 
CMB map combined with the DA's beam window function, plus a noise map computed 
as White Noise x <jq /yNobe, where do is the rms of the noise per number of observations and 
iV bs is the corresponding patch containing the number of observations per pixel. Finally we 
apply the corresponding mask patch. 



B. The convergence maps 

We generate the kernel for each DA. In the calculation of the kernel we use the lensed 
CMB power spectrum measured by WMAP, combined with the window beam function and 
the detector noise characteristic of the DA. For the fiducial unlensed power spectrum we use 
the power spectrum computed with CAMB. We plot the m = mode in Fig. [3| All kernels 
have similar size and similar structures within the 1% level, differing only in the smaller 
scale structures (of order < 0.1°) at sub-percent level. For the numerical implementation 
we truncate the kernel at the 1% level, which corresponds to kernel = 1-1° similarly to the 
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Figure 3: Kernel in real space for the four DA's of the WMAP W-band (obs). We plot 
the dominant {m = 0) eigenfunction only. The levels denote the orders of magnitude as fractions of 
the maximum value, delimited by {1, 1CP 1 , 1CP 2 , 10 -3 , 10 -4 }. The kernels in all panels are derived 
from the WMAP TT power spectrum combined with the beam window function of each DA. Top 
left: Wi, Top right: W 2 , Bottom left: W 3 , Bottom right: W A . 



experiments in the previous section. From each patch i of the WMAP map from the W- 
band DA Wj, we compute the corresponding convergence map k l w . For illustration, we show 
three patches of the W\ DA map obtained from projections centered at the galactic coordi- 
nates (long, lat) = {(0, 45°N), (0, 0), (0, 45°S)}, and the corresponding convergence patches 
(Fig. [4]). The contours of the mask are closely reproduced in the convergence patches, which 
demonstrates that the reconstruction is local. 



C. The convergence power spectrum 

To estimate the spurious correlations due to the unlensed CMB, we compute the con- 
vergence map from various realizations of the unlensed patch corresponding to the single 
realization of the observed patch. The average power spectrum of these convergence maps, 
weighted by the inverse of their variance according to Eqn. ([7]), is our estimate of the Gaus- 
sian contribution of the CMB to the convergence power spectrum of the patch. Thus, for 
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Figure 4: Projections of the W\ DA map and KQ85 mask centred at 

(long, lat) = {(0, 45°N), (0, 0), (0, 45°S)} respectively first, second and third row, and the 
resulting convergence patches. Left panels: The masked CMB patches. Central panels: The 
mask patches, with amplitude on the black regions and amplitude 1 elsewhere. Right panels: 
The k — (^) pixcl patches obtained from the corresponding masked CMB patches. 

each patch % and for each DA Wj, we have 

For each patch, we compute the average over the four DA's, which includes four auto- 
correlations and six cross-correlations in both the C£ and the (C^ =0 ) CMB terms, weighted 
by the inverse of the variance of each (C^). w . w t ■ This operation yields (C^) . . We then 

compute the average over all the patches, weighted by the variance of each (C^). , yielding 
the final result C^* . 

We recall that masking a map causes a degradation in the amplitude of the estimated con- 
vergence power spectrum in comparison to the power spectrum estimated from the complete 
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Figure 5: Dispersion of the estimated convergence power spectrum as a function of the 
lensed CMB realization used to define the empirical calibration. The black dashed line is 
the raw estimated convergence power spectrum obtained from the CMB patches after subtracting 
the convergence power spectrum from unlensed CMB realizations. The black filled line is the 
convergence power spectrum obtained after calibrating the black dashed line with the realization 
of Fig. [2j whereas the grey lines are the spectra obtained after calibrating with the various lensed 
CMB realizations. 



map. This amplitude degradation has been quantified for masks of point sources with vary- 
ing density and radius [H] . We correct the convergence power spectrum of each patch by the 
amplitude degradation inferred for the area of the corresponding mask patch, predominantly 
consisting of contiguous large areas. 

Finally we apply the empirical calibration of Fig. [2] (right panel) derived from the simula- 
tions. We remind that the calibration allows to move from the convergence power spectrum 
estimated from one map to the power spectrum estimated from the mean of a distribution of 
maps. In the previous section, this operation introduced a small extra source of uncertainty, 
the term Var s [C^] in Eqn. (15). This term accounts for the error on the mean of the dis- 
tribution and is subdominant as mentioned earlier on. In this section, we are limited to one 
realization that is not necessarily the mean of a distribution of an ensemble of WMAP-like 
realizations. To be conservative, we apply the empirical calibration with an uncertainty 
containing the sample variance per element of the distribution (as opposed to the uncer- 
tainty on the mean), which is propagated into the calibrated power spectrum. This gives 
an important contribution to the error budget, even though the detection is still dominated 
by the estimator noise as in the simulated case. This contribution is illustrated in Fig. [5j 
which shows the dispersion of the estimated convergence power spectrum if corrected by an 
ensemble of empirical calibrations that could be defined by each realization. 

The results, before and after the empirical calibration, are presented in Fig. [6j for both 
the convergence and the deflection field (which is the quantity most commonly used, for eas- 
ier comparison with the results from other detections). After correcting with the empirical 
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Figure 6: Estimated convergence power spectrum from the CMB maps of the WMAP's 

W-band DAs. The grey diamonds are the weighted average of the difference between the power 
spectrum measured from the patches and the power spectrum measured from the simulations for 
each DA and for each patch. The black squares are the same quantity multiplied by the empirical 
calibration. The dotted line is the theoretical convergence power spectrum computed with CAMB. 
The dash-dotted line is the inverse-sum of the variance of the estimator Varf/c^] for each DA. Left 
panel: The convergence power spectrum. Right panel: The deflection power spectrum. 
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calibration of Fig. [2] (right panel), we recover a lensing signal which agrees with the theoret- 
ical prediction. In terms of the lensing amplitude parameter, after applying the empirical 
calibration, including the calibration dispersion in the error budget, we find a conservative 
value of Ai = 0.98 ± 2.50 for the case of full combination of DA's. 



D. Null test: Measurement of the convergence from noise maps 

We proceed to apply the estimator to noise patches synthesised for each DA from the 
corresponding o"o and iV b s . These are the same maps that we have added to the simulated 
unlensed CMB maps. For each DA, we average the convergence power spectrum over the 
different patches and present the results in Fig. [7j 

The results are consistent with the absence of signal, showing that the estimator is free of 
systematics. We notice that this result did not require a calibration, highlighting the lensing 
origin of the spurious correlations, which are not present in this test on noise maps. 



VI. CONCLUSIONS 

In this manuscript we present a direct detection of the weak lensing from the WMAP 
7-year data. We used the real-space variant of the optimal quadratic estimator originally 
defined in Fourier space. 
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Figure 7: Estimated convergence power spectrum from noise maps of the WMAP's 

W-band DAs. The symbols are the convergence power spectrum averaged over the different 
patches for each DA. Diamond: W\. Cross: Wi- Plus: W3. Triangle: W4. The squares are the 
estimated convergence power spectrum from the CMB maps (after the empirical calibration), 
which we include for comparison. The dotted line is the theoretical convergence power spectrum 
computed with CAMB. 



First we applied the estimator to maps synthesised so as to simulate the WMAP W-band. 
The conventional estimator of the Gaussian bias failed to recover the input convergence 
power spectrum, in contrast with its success in Planck simulations. This failure signaled 
the presence of a non-Gaussian bias. A variation of this estimator was conceived which 
successfully accounted for the additional bias and recovered the input convergence power 
spectrum, with a lensing parameter of Al = 1.55 ± 0.93. The detection error is dominated 
by the variance of the estimator which depends on the experimental resolution and detector 
noise. By comparing the results of the two methods, we derived an empirical calibration 
which estimates the additional bias as a function of I. 

We proceeded to apply the estimator to maps obtained from HEALPix cartesian pro- 
jections of the full-sky maps produced by the WMAP W-band. We used the conventional 
estimator of the Gaussian bias, including the correction of amplitude degradation due to 
masking, and the empirical calibration function. We achieved a good agreement of the 
detection with the theoretical prediction, as well as with the detection from the simula- 
tions, yielding Al = 0.98 ± 2.50. This detection is also compatible with previous detections. 
However, the large error renders this detection not statistically significant, similarly to the 
detection from the WMAP data using the harmonic implementation of optimal estimator 
[?]. The detection from the WMAP data using a kurtosis estimator has higher significance 
because it aims for the lensing correlations of interest [5]. 

Nonetheless, our detection serves as a proof of the robustness of the real-space estimator 
in dealing with extended masks and overwhelming experimental noise. We expect that a 
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full-sky experiment with higher resolution and lower detector noise, such as Planck [18J, 
will allows us to make a more significant measurement of the convergence power spectrum. 
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